Pair Wave Functions in a Bose Liquid 
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Pair wave functions (PWF) which are eigenfunctions of the reduced density 2-matrix are consid- 
ered for a homogeneous Bose liquid. With the Bogoliubov principle of the correlation weakening it is 
proved that the distribution of the "dissociated" pair states over momenta is exactly the product of 
the single-particle distribution functions. Thus, the "dissociated" pair states are naturally classified 
as condensate-condensate, condensate-supracondensate and supracondensate-supracondensate ones 
provided the Bose-Einstein condensate exists. The condensate-condensate as well as condensate- 
supracondensate PWF are expressed in terms of the averages of products of the creation and de- 
struction Bose operators. This leads to the simple interpretation of the anomalous averages as the 
"scattering parts" of the condensate-condensate and condensate-supracondensate PWF. It is shown 
that in contrast to the Fermi liquid, the appearance of the anomalous averages for the Bose liquid 
does not necessarily mean that there exist bound states of pairs of particles. The PWF in the 
Hartree-Fock-Bogoliubov (HFB) approach are found. Given the density of the condensate is not 
zero, there are no bound pair states in the HFB scheme. The expansion of the pair correlation 
function in the set of PWF is very useful in order to take into account both short-range and long- 
range spatial correlations. Applications (possible and already realized) of the formalism developed 
are discussed. 



PACS numbers: 05.30. Jp, 05.30.-d, 03.75.Fi 

I. INTRODUCTION 

It is well-known that for the homogeneous Fermi liquid 
the appearance of anomalous averages (i/j' (xi)ip^ (X2)) 
implies that there exists a macroscopic number of bound 
pairs of particles with zero momentum (here ij)'{x) and 
■0(x) are the field operators). In the case of the pair 
condensate like this, the properties of the system change 
radically. In particular, this leads to the superfluidity of 
the Fermi liquid (superconductivity in an electron sys- 
tem, superfluidity in 3 He). 

For a homogeneous Bose liquid, the superfluidity is 
usually associated with the macroscopic number of par- 
ticles themselves in the zero-momentum state (the Bose- 
Einstein condensate). In the Bogoliubov theory for a 
weakly imperfect Bose gas, repulsive interaction between 
particles depletes the condensate, while makes it thermo- 
dynamically stable. In this case the anomalous averages 
(1)1 (ri)ffi (T2)) play a crucial role also, where i?t(r) is the 
field operator creating a supracondcnsate particle at a 
point r, see Eq. (34). 

Parallels between the Bogoliubov transformations for 
boson jj] and fermion systems gave birth to pair the- 
ories being special cases of the HFB theory (see, e.g. 
Ref. H and references therein). In the framework of the 
HFB approach, the gap in the single-particle excitation 
energy appears Q that may be interpreted, by analogy 
with the theory of superconductivity, as an evidence for 
the bound states of pairs of particles. Below we show that 
in the presence of a single-particle condensate, there are 
no bound pair states of bosons in the HFB theory. Thus, 
it has been verified in terms of PWF that the gap is un- 
physical, which has been stressed in many papers BB. 



The concept of wave functions for a group of m parti- 
cles in medium can be rigorously introduced with eigen- 
functions of the m-th-order reduced density matrix, or 
the m-matrix. Indeed, the system of m particles is a 
subsystem of that of N particles. So, its state is not 
pure even in the situation when the system as a whole 
has a wave function. In general a subsystem is specified 
by the density matrix (see, e.g. Ref. p]). The m-matrix 
is of use when we have a noncoherent superposition of 
the m-particle wave functions. Lowdin and Shull used 
wave functions like these for one and two particles to de- 
scribe electron states in atoms, and called them the "nat- 
ural orbitals" and "natural geminals", respectively (see 
Ref. |Q] ) . The concept of PWF for fermions has been used 
by Shafroth and his coworkers in papers on the theory of 
superconductivity (for a review, see Ref. Q). In what 
follows, we shall use the term PWF for the wave func- 
tions of particle pairs, using the notation of Bogoliubov 
who made the most clear presentation of the concept of 
PWF for fermions ||. 

This paper concerns PWF for the Bose liquid below 
the temperature of the condensation, T c . Analysis is car- 
ried out in general form as far as it is possible, with a 
special emphasis on the bound pair states and that part 
of PWF which is responsible for the correlations of par- 
ticles at short distances. Consideration of these corre- 
lations is of special interest when investigating concrete 
boson systems for which the potential of two-particle in- 
teraction $(r) is, as a rule, strongly singular, or, in other 
words, of the hard-core type (i.e. potential goes to in- 
finity like l/r m (m > 3) as r — > 0, and, consequently, 
there is no Fourier transform of it). For example, the 
well-known Lennard- Jones potential is of this form. In 
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this case standard approximations based on the weak- 
coupling pertrubation theory lead to divergencies. For 
example, the statistical average of the interaction energy 
per particle 



d 3 r<I>(r)g(r) (1) 



(here g(r) is the pair distribution function, n is the den- 
sity of the particles) is infinite in the simplest Hartree 
approximation corresponding to g(r) — 1. A reason for 
the divergency is obvious as the short-range spatial cor- 
relations of two particles are not properly taken into con- 
sideration. In order to escape the mentioned trouble a 
real potential is usually replaced by various effective po- 
tentials for which the Fourier transforms exist (for the 
Bose case see the t- matrix approach fljofl , the pseudopo- 
tential with positive scattering length, or the hard-sphere 
model, [0 and the method by Beliaev This re- 

placement implies that the effective potential contains 
all "information" on the short-range correlations of parti- 
cles. A range of validity of this replacement for the many- 



boson systems is discussed in the textbook 
an inhomogeneous gas, in the recent paper 
discussion at the end of Sec 



IV) 



and, for 
14| (see the 
In the lowest order 
in density n of the Bose gas such effective-interaction 
procedures are reduced to the Bogoliubov model with 
the effective potential [for example, in the pseudopoten- 
tial method the Fourier transform of the "bare" poten- 
tial is replaced by <5(fc) = Ana/m, where a is the scat- 
tering length obtained from the two-body Schrodinger 
equation]. However, this approach is found to be ther- 
modynamically inconsistent |lq| , which can manifest it- 
self in various unphysical results. For example, one can 
obtain no depletion of the Bose condensate within the 
pseudopotential approach, if one tries to apply a self- 
consistent treatment fl6|| . Besides, we are not able to 
obtain the correct value of the interaction energy (Q) by 
direct calculations using g(r) ]l7t]. Thus, it looks inter- 
esting and promising to realize an alternative variational 
scheme that is self-consistent from the very beginning, as 
the short-range correlations are properly taken into ac- 
count for all the terms in the expansion of g(r) (31) in 
the set of PWF. Thus, PWF are "regular channels" which 
can be used to cancel the divergencies in a self-consistent 
manner (see discussion in Sees. £y] and 0). The first 
results derived in Refs. demonstrate that the 

formalism of the PWF can actually be employed as a ba- 
sis for the self-consistent variational treatment for a Bose 
gas with the strongly singular interaction potential. 

General properties of the eigenfunctions and eigenval- 
ues of the correlation function (|^) (differing from the 2- 
matrix only by a norm) arc interesting in itself when an 
eigenvalue becomes a macroscopic quantity (off-diagonal 
long-range order Q ) . With the Bogoliubov principle of 
correlation weakening we shall show that eigenvalues of 
the correlation function (|B|) which belong to continuous 
spectrum ( "dissociated" pair states) are expressed as the 



product of the single-particle momentum distributions. 
Thus, the "dissociated" pair states are naturally classified 
as condensate-condensate, condensate-supracondensate 
and supracondensate-supracondensate ones provided the 
Bose-Einstein condensate exists. The condensate- 
condensate as well as condensate-supracondensate PWF 
are expressed in terms of the averages of products of the 
creation and destruction Bose operators. This leads to 
the simple interpretation of the anomalous averages as 
the "scattering parts" of the condensate-condensate and 



condensate-supracondensate PWF (see Sec. til). 

In the next section, the PWF are introduced for boson 
systems by analogy with Bogoliubov's paper ^| , and the 
continuous spectrum of the correlation function (||) is cal- 
culated. In Sec. Ill we consider the general structure of 
the 2-matrix below T c , the PWF in an explicit form and 
their physical interpretation. In Sec. |v| the PWF are 
evaluated in HFB approach as well as in the Bogoliubov 
model, and some applications of the developed formalism 
are considered. Main results are summarized in Sec. [v|. 



II. THE CONCEPT OF PARTICLE PAIR STATES 
FOR BOSONS 



Let us consider a homogeneous system of N bosons, 
and, for simplicity, let the spin of the particles be equal 
to zero. The Hamiltonian is assumed to be invariant with 
respect to translations and the transformation of parti- 
cle momentum — * Pi. The latter implies that the 
Hamiltonian does not change when the canonical trans- 



formation ak — > a-k, 



is performed. In partic- 



ular, this condition is satisfied for pairwise interactions 
with a potential $(r) = <i>(r). A state of the system is 
determined by a density matrix 



P(ri, ...,r^;ri,...,rjv) = ^ io n *£(ri, ...,r N ) 

n 

x*„(ri,...,r^), (2) 

where *^2 n w n = 1, w n > 0, and \&„ are orthonormal 
system of functions being symmetric with respect to any 
permutations of particles. 
The 2-matrix is defined as 

p 2 (r' 1 ,r' 2 ;r 1 ,r 2 ) = / d 3 r 3 • • • d 3 r N 
Jv 

xp(r'i,r 2 ,r 3 , . . . , rjy; ri, r 2 , r 3 , . . .,r N ), (3) 

and can be expressed in terms of the field operators 
and i/it(r) (see, e.g. Ref. M) 



P2 (r' 1 ,r^;r 1 ,r 2 ) = _ (r^ (^WS^i)) 

(4) 
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Here (• ■ ■) = Tr (• • • pi) stands for the statistical average 
over the state p. The density matrix /? 2 is normalized, 
namely 

Trp 2 = / d 3 rid 3 r 2 p 2 (ri,r 2 ;ri,r 2 ) = 1, 
Jv 

therefore, any matrix element (||) has asymptotic be- 
haviour as l/V 2 , in the thermodynamic limit V — > oo, 
N/V = n = const. Therefore, it is more convenient to 
work with the pair correlation function 



F 2 (ti,T2',t' 1 ,t' 2 ] 



(^(r^t^)^)^)). (5) 



The boundary conditions for F 2 22 follow from the 
principle of the correlation weakening at macroscopical 
separations |J: 



(^(r^ra))^)^)) 



(6) 



when 



ri — r 2 = const, r[ — r' 2 = const, — ri| — > oo; (7) 

(^rO^K^fo)^)} (8) 

when 

ri — r[ = const, r 2 — r' 2 = const, |ri — r 2 | — * oo. (9) 

It should be stressed that limits (Q) and (H) are valid ei- 
ther the Bose-Einstein condensation takes place or not. 
In the first case the anomalous averages in the expression 
@ are not zero p3[ . 

As the expression (^|) is a Hermitian kernel, we can ex- 
pand it in the orthonormal set of its eigenprojectors (in 
the Hilbert space): 

F 2 (r lj r a ;riy 2 )=V;M:(r lj r a )* I/ (riy 2 ) (10) 



i.e. the sum of all N v is the total number of pairs. There- 
fore, the non-negative quantities N v can be interpreted 
as the mean number of the pairs in the state v, any pair 
being doubly taken. The ratio w u = N U /[N(N—1)] is the 
probability of observing a particle pair in the pure state 
with the wave function \&„(ri, r 2 ). Here, as one might 
expect, J2 V w u = 1. 

From the definition (||) it follows that 

F 2 (r 1 ,r 2 ;ri,r 2 )= J F 2 (r 2 ,r 1 ;ri ) r 2 ) = F 2 (r 1) r 2 ;r 2 ,r' 1 ) 

hence, \l/ 1/ (ri,r 2 ) = ^l/ i/ (r 2 ,ri), i.e. the PWF for bosons 
are symmetric with respect to permutations. 

For equilibrium states, the momentum of centre of 
mass of a pair of particles q is a quantum number corre- 
sponding to PWF, so the index v can be represented as 
v — (u>,q), where ui stands for other indices. PWF can 
be written in the following form 



*„(ri,r 2 ) = VVq( r i - r z) 



exp[?q(ri + r 2 )/2] 



V 



(12) 



Let us demonstrate this. Consider an equilibrium state 
described by the Gibbs canonical ensemble. As the 
Hamiltonian is invariant under translations, the total mo- 
mentum P is conserved. So, we have 

*n( r i +a, . . . ,rjv + a) = * n (ri, . . . , r N ) exp[£Pa], (13) 

Without loss of generality we put P = 0, and introduce 
new variables 



R=(ri+r 2 )/2, r = r x - r 2 , 



(14) 



and similarly for r[ and r 2 . Rewrite (||) with the help of 
(|l3| ) setting a = R and a = R', respectively. Substitut- 
ing rewritten Eq. (|2|) into (Q) and performing the Fourier 
transformations with respect to r' 3 — R' , . . . , r' N — R', r3 — 
R, . . . , r^r — R we obtain: 



patten, r 2 ) = I^ e ^ R '- R ) Pq (r',r), 



(15) 



where 



/ d 3 nd 3 r 2 |^(n,r 2 )| 2 = 1. 
Jv 



(11) 



Eigenfunctions \E , y (ri,r 2 ), which at the same time are 
eigenfunctions of 2-matrix (^) , are called the wave func- 
tions of pairs of particles, or PWF. From (0), ([l0|) and 
dTH) we get 



d 3 nd 3 r 2 F 2 (n , r 2 ; n, r 2 ) = ( A^ 2 - N) 



where p q is a non-negative Hermitian kernel. Expanding 
it in the set of its eigenprojectors ^^(r')^* q (r) we ar- 
rive at Eqs. ( ]To| ) and (12). Analogous proof can be given 
for the Gibbs grand canonical ensemble. 

Then the expression (fLOh can be written in the form 



fMn^rV,) = - r 2 )^, q « - r' 2 ) 



cj,q 

x exp 



V 



[if(ri + r 2 



-r x -r 2 ) . (16) 



For ip Ut q(r), which can be interpreted as the wave func- 
tion of a pair of particles in the centre-of-mass system, 
from ( jll| ) and (|l2|) we get the following normalization 
condition: 
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d 3 r |^, q (r)| 2 = 1. 



(17) 



As the wave function V'w.qfr) is an EF of the kernel 
Pq(r\r) (see Eq. (15)), it can belong either to a dis- 
crete or a continuous spectrum, so we have two possible 
situations Q: 

(1) V'o; q( r ) — > at t — > oo (bound state of a pair of par- 
ticles). 

(2) VVq( r ) — * V / 2cos(pr) at r — > oo ("dissociated" states 
of a pair of particles corresponding to scattering with rel- 
ative momentum p, the factor y/2 stands for the appro- 
priate normalization when p ^ 0). 

In the first case, uj = i is a discrete index numbering 
bound pair states. We denote V*w,q( r ) = Vq,i( r ) an d keep 
the normalization (p~7|) : 



/ d 3 r K,(r)| 2 = 1. 
Jv 



In the second case 



VVq(r) 



w = p and 

¥>P,q( r ) 



V 



Due to Eq. @ we obtain the normalization 
- / d 3 r |^ P , q (r)| 2 = 1. 



(18) 



(19) 



With the variables (|l4| ) the expression (|l6|) can be rewrit- 
ten in the following form: 



F 2 (n, r 2 ; r 1; r 2 ) = F 2 (r, R; r', R') = £ -^Vq, i (r) ¥ > q , i (r / ) exp[zq(R - R)] 



+ £ -^^, q (r)^,q(r')exp[iq(R' - R)]. 



p.q 



(20) 



In the limit V — > oo we should replace the quasidiscrete momentum sums by the corresponding integrals: 
F 2 (r,R;r',R') = ]T f d 3 q Wi^W^M exp[zq(R - R)] 

+ J d 3 pd 3 q W (p,q)^ q (r)^p, q (r')exp[iq(R'-R)]. (21) 



Thus, one can see from Eqs. ©, © that Vwi(g)d 3 g 
is the number of bound pairs in state i with momen- 
tum q in an infinitcsimally small momentum volume d 3 g, 
and V 2 w(p, q)d 3 pd 3 q is the number of "dissociated", or 
scattered, pairs with relative momentum p and centre- 
of-mass momentum q in infinitesimally small momentum 
volumes d 3 p and d 3 q. 

In the centre-of-mass system the replacement p — > — p 
corresponds to the permutation of particles. So, the fol- 
lowing symmetric relations take place: 



iu(p,q) = w(-p,q), 
V P ,q(r) = Vp,q(-r) = ¥>-p,q( r )- 



(22) 
(23) 



Let us now prove, with the principle of the corre- 
lation weakening, that the distribution of the particle 
pairs over the "scattering" states w(p, q) is expressed 
in terms of the single-particle momentum distribution 
w S i„ 9 (k) = n S(k) + n(fc)/(27r) 3 , where 5(k) stands for 
the 3-dimensional ^-function, n(k) = (a^fik) is the single- 
particle occupation number. Suppose that the tempera- 
ture is larger than the temperature of the Bose-Einstein 
condensation, so, hq = and w S i„ g (k) is a regular func- 
tion, i.e. it does not contain 5-function after the thermo- 
dynamic limit. On the one hand, in the limiting situation 
of @ we have the relation (|§|), which can be written as 



/d. f d"^ 

^^«(Pi)exp[ipi(r' x -ri)] J J^ n (P?) exp[ip 2 (r 2 - r 2 )] 



d 3 qd 3 p 



3 n(q/2 + p) n(q/2 - p) 
(2tt) 6 



exp[ip(r' — r)] exp[zq(R' — R)], 



(24) 



where, passing to the last equality, we introduced the new variables q = Pi + p 2 and p = (pi — p 2 )/2 and used 
notations (fL4|). On the other hand, when (l9j> is true, we have 



r = r 2 - ri 



oo, r- 



oo, R' — R = const, r' r = const. 
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Hence, in this limit we obtain 

</>q,i( r ) °, ^ P ,q( r ) \/2cos(pr). 

Therefore, it follows from (|l]) that in the limiting case 
(||) we have 

F2(rx,r 2 ;r' 1 ,r' 2 ) — > / d 3 qd 3 p w(p, q)2 cos(pr) cos(pr') 

xexp[iq(R'-R)]. (25) 

Further, the Riemann's theorem |^6| used while integrat- 
ing over p and relation (^2|) allow us to rewrite (25) as 

i r 2(ri,r 2 ;ri,r' 2 ) -> / d 3 gd 3 pw(p,q) exp [ip(r' - r)] 

x exp [iq(R' - R)] . (26) 

The right-hand side of Eq. (|||) is equal to that of (|2^) at 
all the values of space variables r = r'— r and R = R'— R. 
Hence, we derive the following equality: 



tu(p,q) = 



n(q/2 + p)n(q/2 - p) 
(2tt) 6 



(27) 



Below the temperature of the Bose-Einstein condensa- 
tion a macroscopic number of particles Nq occupies the 
zero- momentum state: no = {a\ao) /V — Nq/V — const, 
when V — ► oo, and we have in the thermodynamic limit 

(ft(ri)M)) = J d 3 Pl (n 6( Pl ) + 



x exp [ipi(ri - ri)] . 

By analogy with the above proof one can readily be con- 
vinced that the distribution w(p, q) = u>( Pl ~ P2 , pi +P2) 
contains the following parts: 



1) ng*(pi)*(p2), 

2) n <5(pi)— — r +n 5(p 2 )- 



3) 



(2tt)3 
n(pi) n(p 2 ) 
(2tt) 3 (2tt) 3 ' 



(2tt) 3 ' 



(28) 
(29) 

(30) 



The expression ( p8[ ) corresponds to condensate pair 
states, or condensate-condensate pairs (both particle of 
a couple are in the condensate); the terms (29), to 
condensate-supracondensate pairs; and the last term 
(|30|), to supracondensate-supracondensate ones. Thus, 
owing to the expressions (^8|)-(^0|) we obtain the natural 
classification of the "dissociated" pair states. 

We note that the total number of bound pairs Nf, ~ 
V J d 3 qiVi(q) is asymptotically proportional to V 1 , 
while the total number of "dissociated" pairs Nd ~ 
V 2 Jd 3 pd 3 q w(p, q) = N 2 is asymptotically propor- 
tional to V 2 . From an intuitive point of view it is obvious, 
as a given particle in a liquid can form bound states with 



M particles (M — const as V — ► 00), while it forms "dis- 
sociated" states with the other N — 1 — M particles. So, 
we arrive at the following asymptotic equation: 



N b + N d = N(N - 1) ~ N d 



N 2 



Diagonal elements of F 2 (r 1 ,r 2 ;r' l ,T 2 ) depend on squ- 
ared absolute values of the PWF and determine the 
behaviour of the static pair distribution function g(r), 
which can be directly observed in scattering experiments: 



l(T) 



V 2 



N(N — 1) 



F 2 (r,R;r,R) 



After the thermodynamic limit this expression takes a 
form 

SW = ^(E / d 3 qwM\^Ar)\ 2 

^ i 

+ yd 3 pd 3 (?W (p,q)|^ p , q (r)| 2 ^. (31) 

One can easily be convinced that g(r) — ► 1 as r — > 00. 
The behaviour of g(r) for r — > is determined by the 
PWF at small distances. Because of repulsion, the prob- 
ability of finding three particles in a very small volume 
is much lower than that of finding two particles in the 
same volume. Therefore, the short-range behaviour of 
two particles is determined by the ordinary two-body 
Schrodinger equation with the pairwise interaction <I>(r). 
Thus, the PWF for short distances between two particles 
is proportional to the usual wave function (fp\r) of the 
two-body problem P7||: 



^ p , q (r)^C<)(r), 



(32) 



when r — > 0. This implies that for a singular potential 
the PWF goes to zero for r — * 0, and it follows that 
g(r) — > as r — > 0. 



III. STRUCTURE OF PAIR CORRELATION 
FUNCTION AND PAIR WAVE FUNCTIONS FOR 
BOSE LIQUID WITH CONDENSATE 

In order to examine the behaviour of the pair correla- 
tion function (|5|) with a condensate we will employ the 
procedure proposed first by Bogoliubov (!]] and lately jus- 
tified with the help of the principle of correlation weaken- 
ing M : in the calculation of averages p8| of any number 
of the field operators we can substitute the C-number 
V No for the operators clq and a J involved in $(r) and 
Here Nq is the number of particles in the con- 
densate. It should be stressed that, after Bogoliubov, 



5 



this substitution gives asymptotically exact values of the fact that the main contribution to the value of the inte- 
averages in the limit V — > oo and does not imply any gral comes from infinitely large r. Thus, we can put 



approximations. Indeed, we have for any average that 
involves the operator clq/\V : 



lim (— ;= ■ ■ ■ ip(r') 



4>{r) = Vno + i?(r), 
ft(r) = y/nS+P(r), 



(33) 



where 



= lim — 

y^oo V 



lim - / d 3 r(^(r)---^(r')---) 



d 3 r(^(r))(..>(r')---)=V^<---^(r')---). 



(34) 



p#0 



Here we use the definition (ip(r)) = (ip^r)) = 



(no = Nq/V is the density of the condensate) and the ($*(r)) = we obtain 



are supracondcnsatc field operators. Substituting the ex- 
pression d33| ) into (|^) and using the relation (f?(r)) = 



F 2 (r 1) r 2 ;ri,4)=F 2 (1) (r 1 ,r 2 ;r' 1 ,ry + ^(n,^;^,^) + ^(ri,^;^,^) + ^(n,^;^,^), (35) 
where we introduce the notations 

nl + no0Hn)&(*2)) + "o(i?(r 2 )i?(ri)) ; (36) 



A3), 



^ 2) ( ri ,r 2 ;r;,r 2 ) 
^(n.rajriya) 



^((^(n)^^)) + (^( ri )^(r 2 )) + (^(r 2 )tf(r 2 )} + (tft(r 2 )tf(ri)>j , 
V^((^(r 2 )^(r 2 )^(r' 1 )) + (tft^)^)^)) 

+ (^(n)^)^)) + (tffrj&fatfft)) 
(^(n)^(r 2 )iJ(i4)tf(ri)). 



(37) 

(38) 
(39) 

The expressions (|B6|)-(|39|) can be expanded in systems of projectors (in the Hilbert space). For (|36[)-(|38|) the corre- 
sponding projectors can be expressed in an explicit form in terms of averages of the creation and destruction Bose 
operators a p and ap. For ( |36"| ) we have 

F^n^rV,) = nfo*(n - r 2 Mr[ - r 2 ) - (^Jtftfo))^)^)), 

where 

ip(r) = l + ip(r), 

4>{r) = V>(n - r 2 ) = V*(ri - r 2 ) = (tf (r^tf (r 2 ))/n = ^— / d' 1 /.- < •(/,')< ' kl '. 
V'(fc) = (a k a_ k )/"o- 
In a similar manner we get 



(2tt) 3 



Here 



F 2 (2) (r,R;r',R') = J d 3 pd 3 g w (1) (p, q)\/2cos(pr)V2cos(pr') exp[£q(R' - R)]. 



^W(p,q) = 2n <5(q/2 - p) " (< ^ 2 + P) = 2n 5(q/2 - p)- n(<?) 



(2tt) 3 



(2tt) 3 



(40) 
(41) 
(42) 

(43) 
(44) 



is the momentum distribution for condensate-supracondensate pairs corresponding to the expression (G3) . 



(3) 

Performing the Fourier transformations, we reduce F 2 to the following form 

Fi 3 ) (r ,R; r ',R) = /dV 3 , ^(p,q) (^)V2cos(p r ) + ^ (r) ^cos(p r 0) X exp[iq(R' 



R) 



(45) 



G 



where 



Mr)= 



(2tt) 3 



d 3 k 4> p (k)e 



kr 



^ p (k)= Vp(-k) = ^_p(k) = v^Vo 



1 («2p a P+ka P -k) 

\/2n n(2p) 



(46) 
(47) 



The second equation in ( |47| ) is a consequence of the invariance of averages with respect to the transformation — > 
a_k, & k — > a_ k , as the Hamiltonian does not change when the transformation is performed. The factor yA^o 

compensates the decrease of the average value (a 2p a p+ ka P -k) ~ V^ 1 / 2 , so, the quantity i/> P (k) is of order 1 in the 
limit V — > oo. 



From (46) and ( [47|) we obtain the symmetry properties 

ip P (r) = ^p(-r) = V'-pW- 
Summing the expressions ([43|) and (B^) we get 



F 2 (2) (r,R;r',R')+F 2 (3) (r,R;r',R') = / d 3 pd 3 q itf^p, q)y>*(r)y> p (r') exp[iq(R' - R)] 



where 



- / d 3 < Z 2n ^< /2 (r^ q / 2 (r')exp[ i q(R'-R)], 



¥> P ( r ) = ^-p( r ) = M~ r ) = V2cos(pr) + ^ p (r), 



(48) 



in accordance with (f23|). 

Now it is not difficult to express (p3) in a form similar to (Eh) 



/Vr.R:r'.R') = (D,;(r')+ / d 3 pd 3 <Z 2n 5(q/2 - p ) Jl^ p (r) Vp (r') exp[iq(R' - R)] + F 2 (r, R; r', R'), (49) 



where 



F 2 (r, R; r', R) = F 2 (r 1; r 2 ; r' 1: r' 2 ) = F 2 (4) (n, r 2 ; r' 1; r' 2 ) - (^(r^t^))^)^)) 



- J d 3 q 2n ^V; /2 (r)V q /2(r') exp[iq(R' - R)]. 



(50) 



As it is proved in Sec. [ll], the eig envalues of 
F 2 (r,R;r',R) are given by Eqs. So, 
one can conclude that the projectors tp*(r)<p(r') and 
^ q / 2 ( r ) ex P[~* f l^'] l y 9 q/2( r ') exp[iqR'] are nothing else 
but the eigenprojectors for the pair correlation func- 
tion, and, hence, cp(r) and <^q/ 2 ( r ) exp[jqR] are the PWF 
for the Bose system below T c . The first term in ( [49| ) 
corresponds to pairs in the condensate; the second, to 
condensate-supracondensate pairs; the last one, F 2 , to 
pairs in the supracondensate, and, maybe, to bound 
pairs. The PWF (p(r) and <£> p (r) given by expressions 
(|40|)-(42) and ( |h3)-(|48|) respectively are normalized in 
accordance with (19[). The condensate PWF differs es- 
sentially from 1 at short distances and goes to 1 when 
distances between particles in a pair increase. The sit- 
uation is similar to the scattering of two bare particles 
with zero momenta at infinity. The anomalous averages 



(■(? t (r 1 )^t( r2 )) in PWF @ have essentially different 
physical interpretation in comparison with the fcrmion 
systems: they determine the "scattering" part of the con- 
densate PWF, while for fermion systems anomalous av- 
erages (^(ri)^^)) themselves are PWF for the bound 
state of pairs of particles with zero momentum . 

In the condensate-supracondensate PWF (fl8|) we drop 
the index corresponding to the centre-of-mass momen- 
tum as q = ±2p. This is a consequence of the 5-function 
being present in the pair momentum distribution (Q). 
For the same reason in the condensate-condensate PWF 
d40| ) the indices are omitted at all. 

The functions ip{r) and ip p (r) defined by Eqs. ( (4l| ) 
and (|4^) are responsible for particle short-distance cor- 
relations, therefore, owing to ( |3^ ) for strongly singular 
potentials we have 



V>(r = 0) = -1, ^ P (r = 0) = -V2, 



(51) 
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(see the discussion at the end of Sec. [Hj). 

It should be emphasized that the pair correlation func- 
tion for supracondensate particles is F 2 rather than F 2 
defined by Eq. fl3!3|). Since the kernel flls| ) is non- 
negative, the kernel (^) is also non-negative. Conse- 
quently, the last two terms in the expression for F 2 are re- 
duced to zero by the terms in . These terms look like 
those corresponding to bound states of pairs. If the kernel 
F^ written in form (till) contains other terms associated 



with bound states, then the latter are also contained in 
F 2 and, consequently, in the final expression (49). The 



distribution over supracondensate-supracondensate pair 
states is given by the expression (|o|), and thus, we ar- 
rive at the final formula for F2 : 

F 2 (r,R;r',R') 

= Y,j d 3 9 ^(g)^(r)^ q ,,;(r')exp[ iq (R'-R)] 

i 

,3 ,3 n(q/2 + p) rt(q/2 - p) 

dpdq ^f — wr~ 

x ^, q (rV P ,q(r') exp[iq(R' - R)], (52) 



where the symmetry properties (23) are fulfilled. Thus, 
in order to answer the question about real bound states of 
pairs of particles, one must carry out a detailed analysis 
taking into account all groups of the terms (|36|)-(|3"9|). We 
emphasize that we do not specify the size of the bound 
pairs [the characteristic range of the PWF y q ,i(r)]. If 
it is much more than the mean distance between parti- 
cles, then, following Bogoliubov ||, one may call these 
pairs "quasi-molecules" pjjj . If the radius of the bound 
particle couples is of the order of the mean distance be- 
tween particles or, even, less then one may speak about 
ordinary molecules. Note that "quasi-molecules" are not 
stable for a Bose system within the mean-field variational 
treatment |3(]] which is equivalent to the HFB approach. 
It is a difficult question whether it is the case beyond the 
mean-field approximation, so, a priori we cannot ignore 
existence of "quasi-molecules". The bound pair states 



can totally be a result of the collective effects, as in the 
case of Fermi systems J3lj . We note that such situation 
can even be realized in the case of three particles |32| . 

It is not difficult to see with the help of the Riemann 
theorem p6| that the boundary condition (Jfy is satisfied 
for ([H, f>|) in the limit (@). 

The expressions (|3^)-(|39[) can be interpreted in the 
conventional way, in terms of scattering processes for 
pairs of particles (pi,P2) — > (p'^p^), where pi, P2, 
p[ and p 2 are the particle momenta. For this pur- 
pose one should calculate the averages in the momen- 
tum representation, i.e. make Fourier transformation. 
This procedure corresponds to expansion of the PWF 
in plane waves. Thus, the first term in (|3^) is related 
to Pi = P2 = Pi = P2 = 0; the second one, to 
Pi, P2 ^ and p[ — p' 2 = 0; and so on. The terms 
in (Bq) correspond to processes when one of the four 
momenta is zero. These processes become dominant in 
the vicinity of T c Q, when Nq/N <C I. In this case 

(^) 2 < ^ < W^f, so one can expect that F 2 (3) is im- 
portant in the expression (^5|) . Taking into account that 
the thermodynamic potential is tightly connected with 
the correlation function, one can assume that the terms 
in (|38|) are important to obtain correct thermodynamics 
at those temperatures. Recall that these are precisely the 
terms being responsible for the short-range correlations 
for the condensate-supracondensate pairs. 



IV. PAIR WAVE FUNCTIONS IN 
HARTREE-FOCK-BOGOLIUBOV AND 
BOGOLIUBOV MODELS 

Let us consider the HFB approach from the point of 
view of PWF. In the HFB method (see the detailed dis- 
cussion in Ref. 0), a Hamiltonian is approximated by 
a quadratic form of the Bose operators a p and a p con- 
nected with initial operators a p and a p by the canonical 
Bogoliubov transformations. Consequently, all averages 
can be calculated with the Wick-Bloch-De Dominicis the- 
orem Thus, we have 



F 2 (3) (r 1 ,r 2 ;r' 1 ,r 2 ) = 0, 



(53) 
(54) 



The correlation function (p5|) can easily be rewritten in the following form 



F 2 (r,R;r',R) = nfo*(r)<p(r') + d 3 pd 3 q 2n S(q/2-p) 



n(q) n(q/2 + p) n(q/2 - p) 



(2tt) 3 (2tt) 3 
\/2cos(pr)\/2cos(pr') exp[iq(R' — R)], 



(2tt) 2 



(55) 



where ip(r) is related to the averages via (fl0|). In so doing, n(p) and ijj(k) are expressed in terms of Hamiltonian 
parameters in a self-consistent manner. From (|55|) one can see that there are no bound states of particles in the 
framework of the HFB approach. We note that the expression (|55|) corresponds to the ideal Bose gas when (p(r) = 1 
and n(p) is the Bose-Einstein distribution with a zero chemical potential. One can see that in this case all PWF are 
symmetrized plane waves, as it is to be expected. 
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In the HFB method, only the condensate PWF differs from the plane wave. Since PWF are intimately related to 
the excitation spectrum, one can assume that it is a rough approximation, which has led to the unphysical gap in the 
single-particle excitation energy. As Eq. ( |53| ) is fulfilled for all temperatures, then the HFB model cannot correctly 
describe the thermodynamic behaviour of the system in the vicinity of T c (T < T c ). It should be pointed out that the 
expression (H) is valid provided n ^ 0. If n = 0, but (t?t( ri )tft( r2 )) -L q, as in the theory of Valatin and Butler |35| 

(see the careful analysis by Nozieres and Saint James |3(]]), then Frp = — — and F^ is determined by 
(|54|). In this case the anomalous averages really correspond to the bound pairs of particles with q = 0: 

F 2 (r,R;r',R') = W *(r)^') + J d 3 P d\ w( ^+ P) n(< *^~ p) y/2 cos(pr) ^2 cos(pr') e X p[iq(R' - R)], 



where pa is the density of the bound pairs 

P o = [ d 3 n IWn)^))! 2 , 



and <p(ri — r 2 ) = (^(ri)^^)) A/po is PWF normalized 
in accordance with ( |l8| ) p6[ . 

The Bogoliubov theory of weakly imperfect Bose gas is 
a particular case of the HFB theory when the condensate 
depletion is small: (n — n^)/n <C 1. It can be obtained 
directly from (|35| ) if we neglect terms of the third and 



forth orders in d and ir 37 



^ 2 (r,R; r',R') 

= nl<p*(rMr') + J d 3 pd 3 q 2n <5(q/2 - p)^L 
x\/2cos(pr)\/2cos(pr')exp[iq(R' - R)]. (56) 



Comparing the exact expressions (J49j) and (|52|) with 
Eq. (56) one can reveal that the Bogoliubov model 
involves the following approximations. First, the 
supracondensate-supracondensate part is omitted, which 
is completely justified as the condensate depletion is 
small. Second, the condensate-supracondcnsate PWF are 
symmetrized plane waves h^p(r) = in Eq. (ft8|)l, hence, 
the boundary conditions (p^) for a hard-core potential 
cannot be satisfied by any choice of possible parameters. 
In particular, the change of the "bare" potential by the 
effective one leaves this property of the model unaltered. 
The approximation (fi p (r) — V2cos(pr) implies that the 
Bogoliubov model is nothing else but the ideal gas ap- 
proximation for the condensate-supracondensate PWF, 
while the scattering part of the condensate-condensate 
PWF © is not zero but small and can be evalu- 
ated within the Born approximation for the in-medium 
PWF [|l5 18 3£|. Therefore, the Bogoliubov model does 
not properly take into account the short-range correla- 
tions of bosons. So, in the framework of this model 
we obviously obtain the divergency when evaluating the 
mean energy for the "bare" hard-core potential due to the 
infinite contribution of the condensate-supracondensate 
"channel" involved in the representation ( pl| ) to the 
mean interaction energy per particle (]l|). The hard- 
sphere model can be treated as the Bogoliubov model 
with the replacement $(r) — » <5(r)4-7ra/m. Even af- 
ter this substitution we face another divergency. The 



point is that the Bogoliubov model involves, in the im- 
plicit form, the next-to-Born term for the scattering am- 
plitude: at = -b, b = TO/[47r(27r) 3 ]/ d 3 fc$ 2 (fc)/(2T fe ), 
where = fc 2 /(2m). The hard-sphere replacement im- 
plies that <&(fc) = Ana/m — const, which leads to the 
divergency associated with the parameter b, and, hence, 
to divergent the mean energy. A cure for this difficulty 
is to replace the Bogoliubov expression (56) by the exact 
formula (|49| ) in which the term F2 can be neglected due 
to a small condensate depletion: 



F 2 (r,R;r',R')=n 2 <p*(r)tp(r') 

+2n |^3n( g )^ /2 (r)^ q/2 (r')exp[iq(R'-R)]. (57) 

In so doing the PWF <p(r) and <^ p (r) should be deter- 
mined in a self-consistent manner which provides the 
boundary conditions ([32]) [and, hence, (f5j|)]. Along this 
line with the help of the variational procedure we can 
derive the following non-linear equation: 



U(k) = $(fc) 



d 3 q $(|k - q\)U(q) 
( 27F ) 3 • If* + 2nf q U(q) 



(58) 



valid at zero temperature and small densities |T^ , [l8]| . 
Here for the function T q we can put, within a good accu- 
racy, T q ~ T q = q 2 1 (2m), <I>(fc) is the Fourier transform of 
the potential $(r), and U(k) = J d 3 r ip(r)3>(r) exp(— ikr) 
can be treated as a scattering amplitude in a medium. 
In the limit n — ► Eq. ( |58| ) is reduced to the or- 
dinary Lippmann-Schwinger equation for the scatter- 
ing amplitude which corresponds to the zero scatter- 
ing momentum. Note that Eq. (|58| ) looks like an 
equation for the many-body i-matrix r(k, k', q; lj) at 
uj = (see, e.g. the recent review Sec. 4). How- 
ever, there exist essential differences. First, in general, 
the equation for the t-matrix is frequency dependent in 
contrast to Eq. (|5^). Second, T(k, k',q;tj = 0) can 
not be associated with in-medium scattering amplitude 
^k', q (k) = f d 3 r (ik'.qfO^f) exp(— ikr). In particular, 
the condensate-supracondensate PWF (f46])-(f48|) is char- 
acterized by the index k' = ±q/2 and expressed via 
"triple" averages (^) which are completely neglected 
in the many-body i-matrix approximation as well as in 
the Bogoliubov, HFB and Popov ones |3S|]. Thus, one 
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cannot associate the PWF </Jk',q( r ) with the function 
Xk',q:w=o( r ) whose Fourier transform is usually defined 
as rfk,k',q;w) = 1/(2tt) 3 / d 3 p $(p)xk-. qiW (k - p) [see 
Ref. ||],Eq. (A.23)]. The latter is often called "effective 
wave function in a medium". Third, Eq. (58) is obtained 
by means of variational method beyond a mean-field ap- 
proximation Jl5|,[l8| , while an equation for i-matrix is usu- 
ally derived by summing a certain set of the diagrams. 
For this reason the relation T q ~ T q = q 2 /(2m) is an 



approximation for Eq. (|38|) while the i-matrix formalism 
deals exactly with T q . Thus, the relation between the t- 
matrix and PWF formalism is rather subtle and beyond 
the scope of this paper. We only note that our defini- 
tion of the PWF as the eigenfunctions of the 2-matrix is 
consistent with that of quantum mechanics. Namely, we 
operate with the functions which determine the average 
energy of the interaction per particle ([!]) by means of the 
representations fl3l] ) in accordance with the formalism of 
quantum statistical mechanics. While the physical sence 
of the functions Xk'.q;u=o( r ) is rather unclear. 

At sufficiently large momentum k the main contribu- 



tion in the integral in the right-hand side of Eq. J58j) 
comes from the large momenta q. As lim g _,oo U(q)/T q 
0. for k — ► oo Eq. (|58| ) is also reduced to the two- 
body Lippmann- Sch wing er equation, which leads to the 
boundary condition (p2[). However, it follows from Eq. 
(|58|) that at finite density n tp(k) oc 1/k at small k, and, 



hence, ip(r) oc 1/r at r — > oo, in contrast to the two- 
body problem which implies ip{r) oc 1/r. Thus, peculiar 
overscreening takes place for the condensate-condensate 
PWF. Details concerning Eq. (|58|) can be found in 
Refs. HHII. 

With the help of the formalism of the PWF one can 
obtain the exact relationship between the chemical po- 
tential and the PWF J40|) and ( [l8| ) in the presence of the 
Bose-Einstein condensate. This relationship is of special 
interest, for it allows us to see in an explicit form how the 
renormalization of an interatomic interaction takes place 
as applied to the Gross-Pitaevskii equation, at least in 
the homogeneous case. Let us consider at first an inho- 
mogeneous Bose system with the Hamiltonian 



H = I d 3 n 



^ (r ' t] (" 2^ + Vext[r) ) ^ (r ' t)+ d3rdV $(|r " r ' l) ^ t (r ' ( r '< ^ r '' *)' 



(59) 



where the potential V ex t(r) corresponds to the external forces, and 5>(r) is the "bare" interaction potential. Note that 
the Bose field operators used in Sees. || and [II are related to the zero time: ip(r) = ip(r, t = 0), ^(r) = ip^(r, t = 0). 
Using the Heisenberg equation of motion with the Hamiltonian ( |59| ) for the Bose field operator $(r, f) and taking 
average values we get 



■ 00(r,t) 

dt 



2m 



+ V«rt(r)U(r,t)+ / dV$(|r-r'|)(^(r',t)^(r',^(r,t)), 



(60) 



where (f>(r,t) = (tp(r,t)) [and, respectively, cf>*(r,t) = (^>t(r,t))]. On the other hand, one can make use of the 
Bogoliubov's procedure of "extracting" the c-number parts of the Bose field operators in the Hamiltonian (|5^) (see 
discussion in Sec. p|): 



4>(r,t) =<Kr,i)+tf(r,*), 
$ j {r,t) = 0*(r,i) +tf t (r,f). 



(61) 



Thus, in the stationary equilibrium case the grand canonical potential Q depends on the c-number complex parameter 
(/>(r) that can be considered as the order parameter. So, the equilibrium value of the grand canonical potential 
Q = fi(T, V, fx; {(f>(r)y, {</>*(r)}) corresponds to a minimum with respect to the order parameter (the parameter \x 
stands for the chemical potential). Using the well-known expression for an infinitesimal change of the potential 
Sfl = (S(H — /.iN)) one can obtain from the condition 5£l/5(f>(r) = 5£l/5<fi*(r) = 0: 



dV $(|r - r'lXV'Hr'Mr'Mr)} 



(62) 



It is not difficult to see that Eq. fl62j ) is nothing but the stationary form of Eq. ( J60| ) corresponding to the solution 
(f)(r,t) = 0(r) exp(— ifit) because in this case we have for the time-dependent "triple" averages 



(^( r2) t)^( r ^^( r ' 1 ,t))= eX pH M t)(^t( r2 )^(4)^( r ' 1 )). 
Indeed, the principle of the correlation weakening requires that at | r 1 1 — > oo 

(^(rO^ra)^)^)) = (^(n,t)^(r a ,t)^(i^,t)^(ri,t)) 



(63) 



(^ t (ri 1 *))<^(r 2 ,*)^(^,i)^(ri 1 t)) = (^t( ri ))^t( r2 )^( 1 /)^( r /)). 
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The relation ( |63| ) follows from the last equation as 
$ + (ri)) = 0*(n) and $t( ri)t )) = 0*(n) exp(£ M t). The 
stationary Gross-Pitaevskii equation can be obtained 
from the exact equation ( |62| ) with the help of the re- 
placement: 

$(r) -» — (Jfr), (64) 
m 

(^t(rOV'(r')^(r)) -> 0*(r , )^(r / )0(r), (65) 

where a is the scattering length which can be defined with 
the help of the "bare" scattering amplitude: 4ira/m = 
f/(o)( ) [[/(°)(fc) is the solution of Eq. © at n = and 
T q = T q = q /(2m)]. However, within that approach it is 
not clear how the renormalization (p4| ) takes place. The 
most transperent way of understanding this replacement 
in the homogeneous case is based on the PWF formal- 
ism. In this situation V ext {r) = 0, <f> = ^/tiq cxp(zx), 
and we can put the phase x equal to zero. So, for the 
homogeneous system the exact equation (p2) reads 



1 



dV *(|r - r'|)(V» t (r')^(r')V(r)). (66) 



This equation has been obtained in the momen- 
tum representation in the Bogoliubov's paper || (see 
also Ref. po[ ). Using the specific expressions for 
the scattering parts of the condensate-condensate and 
supracondensate-condensate PWF given by Eqs. ( (42| ) 
and ([47|) , one can represent Eq. (Ejff) in the following form: 



H = n o U(0) + V2 



d 3 q 



n(9)^ q / 2 (q/2), (67) 



here t/ p (k) = J d 3 r </j p (r)$(r) exp(— ikr). Equation 



(|67| ) is the exact representation for the chemical po- 
tential via the condensate-condensate and condensate- 
supracondensate PWF and the occupation numbers n(k). 
It can be used in order to obtain the density expansion 
for fi at n — ■> 0. With the help of the expression for 
n(k) established within the variational scheme mentioned 
above 15 IE 




T k + nU(k) 
Tl + 2nT k U(k) 



- 1 



(68) 



one can derive for the condensate density 

d 3 fc n{k)' 




(69) 



where the substitution k = \t'\/2mn in the integral is 
employed. Besides, at n — > from Eq. (58) one can ob- 
tain the density correction for the in-medium scattering 
amplitude |T§|]| 



17(0) = U (Q) (0)[ 



(70) 



Using Eqs. (69) and ( [70| ) and making the "scaling" sub- 
stitution q = qV2mn in the integral in Eq. ( |67| ) 
arrive at the familiar density expansion 11Q-IT2] 



wc 



V = nU^{Q) ^l + ^vW+---j . (71) 

In so doing, the explicit formula for the condensate- 
supracondensate PWF y p (r) are needless, but only the 
relation lim p ^o fp( r ) — V^pi 1 ") is of use. The factor v2 
appears as due to the normalization condition (19). 



Now it is not difficult to see that in the homogeneous 
case the Gross-Pitaevskii equation corresponds to the 
simplest Hartree approximation with the renormalized 
potential @): 



tf<°>(0). 



(72) 



The range of validity of this approximation is determined 
by the next-to-leading terms in Eqs. (|69|)-([7l|) and read- 
ily expressed as v 



', 3 <C 1. In the paper of Prokakis et 
al. |14j the criterion for validity of the Gross-Pitaevskii 
equation is formulated as 



4:Trh 2 a 



(73) 



where the frequency ui characterizes the harmonic trap 
potential. In the homogeneous case (ui = 0) it cannot be 
fulfilled, while the Gross-Pitaevskii approximation (pi]), 



(|65|) is valid, if the condition 



<C 1 is realized. Thus, 



the criterion ([73]) appears to be refined. 



V. SUMMARY 

The reduced density matrix of the second order (^J) is a 
fundamental characteristic of a many-particle system, its 
eigenfunctions, PWF, being the pure states of two par- 
ticles, selected in an arbitrary way, and its eigenvalues 
being the probabilities of finding a couple in those states. 
Thus, the properties of the 2-matrix are of interest in it- 
self and from the point of view of various applications as 
well. 

With the Bogoliubov principle of correlation weaken- 
ing we express the momentum distribution of "dissoci- 
ated" pair states in terms of the single-particle occupa- 
tion numbers [see the expressions (pq)-(|3b])]. This allows 
us to represent the pair correlation function (||) in the 
form (49), (p4), where the condensate-condensate and 
condensate-supracondensate PWF are explicitly given by 
Eqs. ©-@ and respectively. It should 

be stressed that this representation is exact and do not 
implies any mean-field approximation as well as various 
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model assumptions since it is solely based on the Bogoli- 
ubov's principle of the correlation weakening. The trans- 
parent physical interpretation of the anomalous averages 
( |4l| ) and (Q) as the "scattering parts" of the correspond- 
ing PWF is proposed. The boundary conditions (|5l]) are 
obtained in the case of hard-core, or strongly singular, 
potentials. The condensate-supracondensate PWF is de- 
termined by the "triple" average ( |47j ) . The importance of 
these averages is emphasized in Rcf . |l4| for the effects of 
condensate-exited-state interactions in the evolution of 
condensate mean field, which are in accordance with our 
results. In the textbook jllj] the crucial role of the contri- 
bution of the "triple" averages is discussed in the context 
of "dielectric formalism" approach. Those averages are 
not zero in so called Beliaev-Popov approximation dis- 
cussed recently in fl39| . 

Thus, the anomalous averages (#(ri)$(r2)) are re- 
sponsible for the spatial correlations of two particles in 
the condensate, and do not imply that bound states of 
pairs of particles exist; the terms (|38|) are found to be 
responsible for the spatial correlations of condensate- 
supracondensate pairs of particles and can be expected to 
be important in thermodynamics near the critical tem- 
perature T c . In the HFB approach, the PWF are sym- 
metrized plane waves except for the condensate pairs pro- 
vided the condensate density is not equal to zero. 

The formalism of PWF is useful to construct various 
approximations which take into account both short-range 
and long-range spatial correlations of particles in a self- 
consistent manner. This is due to the fact that all the 
pair wave functions in the expansion (^9|), ((5^) are "reg- 
ular channels" from the point of view of its behaviour 
at small r. Along this line one can obtain the system of 
non-linear intcgro-differential equations for PWF, which 
in the weak-coupling approximation leads to the Bogoli- 
ubov model of the weakly interacting Bose gas [pSf , For 
a dilute Bose gas the non-linear integral equation ( p8| ) 
that takes into account both short- and long-range spatial 
correlations can be obtained with the help of variational 
method beyond a mean- field approach [|l7| [l9|. The ex- 
act relationship between PWF and the chemical poten- 
tial (|6?]) allows us to obtain the density expansion for 
the chemical potential (^Tj) in the strong-coupling regime 
(Bose gas is dilute, but the potential is strongly singular 
as the Lennard- Jones one) . In this regime one can derive 
the expansion for the mean energy per particle in powers 
of the boson density n, the famous results jll],[l2| being 
reproduced without any divergencies which are inherent 
in the pseudopotential approximation Jl7| |lSf| . T hus, this 
expansion can be readily obtained from Eq. (|7l|) with the 
thermodynamic relation fi = d(e(n)n) / dn: 



n 

s = — dn' fi(n') 
n J 



ticle (Q) (as well as the kinetic one) [[TtJ. However, this 
is impossible within the effective-potential approach of 
Refs. [^0|-p^[, which obviously indicates the incorrect be- 
haviour of the pair distribution function g{r) (see the 
detailed discussions in Refs. [ ]l7| , |l9| | ) . 

It is worth noting that the concept of PWF is help- 
ful not only in the case of the Bose systems but also in 
the situation of the Fermi ones. Indeed, employing this 
formalism one can easily prove that below the tempera- 
ture of the superconducting phase transition there always 
exist the bound states of fermion pairs beyond the pair 
condensate M . 

The author would like to thank A. A. Shanenko for 
fruitful discussions. 



Moreover, in the framework of our approach we can ob- 
tain by direct calculation the interaction energy per par- 
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